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An effect of large permanent charge: Decreasing flux to 
zero with increasing transmembrane potential to infinity 

Liwei Zhang* Bob Eisenbergf and Weishi Liu* 


Abstract 

In this work, we examine effects of large permanent charges on ionic flow 
through ion channels based on a quasi-one dimensional Poisson-Nernst-Planck 
model. It turns out large positive permanent charges inhibit the flux of cation 
as expected, but strikingly, as the transmembrane electrochemical potential for 
anion increases in a particular way, the flux of anion decreases. The latter phe¬ 
nomenon was observed experimentally but the cause seemed to be unclear. The 
mechanisms for these phenomena are examined with the help of the profiles of 
the ionic concentrations, electric fields and electrochemical potentials. The un¬ 
derlying reasons for the near zero flux of cation and for the decreasing flux of 
anion are shown to be different over different regions of the permanent charge. 
Our model is oversimplified. More structural detail and more correlations be¬ 
tween ions can and should be included. But the basic finding seems striking and 
important and deserving of further investigation. 


1 Introduction 

Membranes define biological cells. They provide the barrier that separates, and de¬ 
fines the inside of a cell from the rest of the world. Membranes are much more than 
just a barrier. They provide pathways for selected molecules to enter and leave cells. 
The barriers must not be perfect or cells would soon die from lack of energy or drown 
in their waste. Biological cells need energy to survive and that is provided (in almost 
all cases) by substances that must cross the membrane. 

Substances cross membranes through proteins specialized for the task. For a very 
long time ME]) these proteins have been separated into two classes, channels 
and transporters (jiOj), and studied in two traditions, one of electrophysiology ( [T8l 
146]). the other of enzymology ([T9, [36]), although the distinction between the two 
approaches was less than clearcut an])- 

Channels have been viewed fundamentally as charged ‘holes in the wall’ (created 
by the membrane) that could open and close (|35]) but, once open, the channel 
followed simple laws of electro diffusion ([8]). 
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Transporters were viewed as more biological devices, involving conformation changes, 
coupling to energy sources (either ATP hydrolysis or the movement of other ions), 
with quantitative description much more difficult, particularly if the description was 
to be transferrable with parameters that were independent of conditions. 

The enormous literature can be sampled in [T] 12] 13] 136] I4H] . Structural biology has 
shown that transporters and channels have very similar structures m eh m\ su). 
Biophysics has shown that the processes that open and close channels (‘activation’ 
and ‘inactivation’) can be coupled to give properties rather like transporters. Physics 
has shown that the electric fields assumed to be constant in classical electrophysiology 
must depend on the distribution of charge in the channel and surrounding solutions 
([9]) and so must change with experimental conditions and with location. 

The detailed properties of open channels have not been compared with trans¬ 
porters in the modern literature, as far as we know, certainly not using models that 
satisfy the physical requirement that potential profiles (i.e., electric fields) be com¬ 
puted from (and thus be consistent with) all the charges in the system. 

Here we consider a simple model of a permanently open ion channel. (We leave 
gating for later consideration.) Most biologists imagine that if the driving force for 
electrodiffusion is increased-that is to say, if the gradient of electrochemical potential 
across the channel is increased in magnitude-the flux through the channel should 
increase. We show here that is not always the case. Consider a channel with large 
permanent charge and the flux of ions with the opposite sign as the permanent charge 
(called counter ions in the ion exchange literature ([17j) or majority charge carriers 
in the semiconductor literature (pi, -38l S3])- The flux of ions with the opposite sign 
as the permanent charge in a channel can decrease dramatically as the driving force 
increases - we term this phenomenon as the declining phenomenon. More precisely, 
if the concentration of the ion is held fixed on one side of the channel, and the con¬ 
centration decreased on the other (‘trans’) side of the channel, the flux of the counter 
ions can decrease if the permanent charge density is large, as we show here. A de¬ 
pletion zone can form that prevents flow even though the driving force increases to 
large values. It is worthwhile to emphasize that if one increases the transmembrane 
electrochemical potential in a different manner, such as, by increasing the transmem¬ 
brane electric potential or the concentration of the ion at one side of the channel, 
then one does not have the declining phenomenon (see Remark 15.11) . 

The decline of flux with trans concentration has been considered a particular, even 
defining proerties of transporters, involving conformation changes of state and other 
properties of proteins less well defined (physically) than electrodiffusion Q2U 42]). 
Declining flux has been called exchange diffusion (in the early transport literature) 
and self-exchange more recently and is an example of obligatory exchange. Obliga¬ 
tory exchange is a wide spread, nearly universal property of the nearly eight hundred 
transporters known twenty years ago mm) with many more known today ([39]). 
Obligatory exchange is ascribed widely to changes in the structure of transport pro¬ 
teins, to conformation changes in the spatial distribution of the mass of the protein 
(P3D- Obligatory exchange is often thought to be a special property of transporters 
not found in channels. 

The structure of many transporters is now known thanks to the remarkable ad¬ 
vances of cryo-electron microscopy, recognized in the 2018 Nobel Prize. A transporter 
(of one amino acid sequence and thus of a perfectly homogeneous molecular type) ex- 
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ists in different states. Each state is said to have a different conformation meaning, in 
physical language, that the spatial distribution of mass is different in different states, 
and the distributions of the different states form disjoint sets, with no overlap. The 
movement of ions is not directly controlled or driven by the conformation of mass, 
however. Rather the distribution of mass produces a distribution of steric repulsion 
forces, and a spatial distribution of electrical forces (because the mass is associated 
with charge, mostly permanent charge of acid and base groups of the protein, but 
also significant polarization charge, as well). It is the conformation of these forces 
that determines the movement of ions. The spatial distribution of forces contributes 
to the potential of mean force reported in simulations of molecular dynamics. 

This paper shows that channels with one spatial distribution of mass can have 
properties of self-exchange (for majority charge carrier counter ions) if the density 
of permanent charge is large. The spatial distribution of electrical forces can change 
and create a depletion zone that controls ion movement, while the spatial distribution 
of the mass of the protein does not change. The conformation governing current 
flow is the conformation of the electric field - the depletion zone - more than the 
conformation of mass, in the model considered here. The current flow of counter ions 
is much greater than the current flow of co-ions because there are many more counter 
ions than co-ions near the permanent charge. Transporters almost always allow much 
less current flow than channels. 

It should be emphasized that the depletion zone considered here arises from the 
self-consistent solution of the Poisson-Nernst-Planck (PNP) equations of a specific 
model (large permanent charge, counter ion transport) and that parameters are not 
adjusted in any way to create or modify the phenomena. This is in stark contrast 
to calculations of chemical kinetic models that involve many adjustable parameters, 
without clear physical meaning, and equations that do not conserve current ([6]). 

Depletion zones play crucial roles in the behavior of nonlinear semiconductor 
devices although there they usually arise at locations in PN diodes where permanent 
charge (called doping in that literature) changes sign. Depletion zones of the type 
studied here are likely to occur in semiconductors but have received little attention 
because they have less dramatic effects than those in diodes ((34; 3£l [43] ) that follow 
drift diffusion and PNP equations rather like those of open ionic channels mm)- The 
possible role of depletion zones in channel function has been the source of speculation 
and experimental verification ( (lOl 32] ). It is striking that depletion zones can change 
the conformation of the electric field and mimic the obligatory exchange traditionally 
thought to occur only in transporters. Depletion zones can create plastic electric 
fields whose change in shape dominate the flux through a channel of fixed structure. 

Our model is of course oversimplified as are any models, or even simulations 
in apparent atomic detail, of condensed phases. More structural detail and more 
correlations between ions can and should be included. But the basic finding that 
large permanent charge can produce depletion zones and those regions can produce a 
decline of counter ion flux when driving forces increase seems striking and important 
and deserving of further investigation. 

The rest of the paper is organized as follows. In Section [2] we describe the quasi- 
one-dimensional Poisson-Nernst-Planck type model and its dimensionless form for 
ionic flow. Assumptions are specified with two key assumptions: a dimensionless 
parameter e defined in (12.71) is small and a dimensionless parameter Qo > 0 defined 
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in AMD from the permanent charge is large. In Section [3l approximation formulas 
(in small e and v = l/Qo) f° r fluxes are provided, which have a number of non¬ 
trivial consequences. One of the apparent non-trivial consequences is that the leading 
order term J\q of cation flux is zero, independent of the values of transmembrane 
electrochemical potential of the cation. The mechanism of this distinguished effect 
of large (positive) permanent charge is examined in details in Section [4] with the 
help of the internal dynamics (profiles of the electric field, cation concentrations 
and electrochemical potential). It turns out the mechanism for Jio = 0 is different 
over different regions of permanent charge. In Section El a rather counter-intuitive 
declining phenomenon - increasing of anion transmembrane electrochemical potential 
leads to decreasing of anion flux - is shown to be consistent with our analytical result. 
Thus, for the first time (to the best of our knowledge), a mechanism of large permanent 
charge for such a phenomenon is revealed. The mechanism is then examined in details 
again with the help of the internal dynamics. We conclude this paper with a general 
remark in Section El 

2 Classical PNP with large (positive) permanent charge 

2.1 A quasi-one-dimensional PNP model 

Our study is based on a quasi-one-dimensional PNP model for open channels, first 
proposed in [33] and, for a special case, rigorously justified in [28]. For a mixture of 
n ion species, a quasi-one-dimensional PNP model is 

J ^^ ri X)s 0 A (X)§)=-e 0 (±.,C, + Q { X)), 
y ' «=i (2.1) 

w =0 ' - J * = * = 

where X € [ao, &o] is the coordinate along the axis of the channel and baths, A(X) is 
the cross-sectional area of the channel at the location X, eo is the elementary charge, 
£o is the vacuum permittivity, e r (X) is the relative dielectric coefficient, Q(X) is the 
permanent charge density, ks is the Boltzmann constant, T is the absolute tempera¬ 
ture, <!> is the electrical potential, and, for the fcth ion species, Ck is the concentration, 
Zk is the valence, X\.(X) is the diffusion coefficient, //*, is the electrochemical potential, 
and Jk is the flux density. 

Equipped with system m , a meaningful boundary condition for ionic flow 
through ion channels (see, [12] for a reasoning) is, for k = 1, 2, • • • , n, 

<h(ao) = V, Cfc(oo) = Ck > 0; d>(6 0 ) = 0, Ck(b 0 ) = 7 Zk > 0. (2.2) 

Mathematically, we will be interested in solutions of the boundary value problem 
(BVP) (12.11) and (12.21) . An important measurement for properties of ion channels is 
the I- V (current-voltage) relation where the current X depends on the transmembrane 
potential (voltage) V and is given by 


X=Y J ZsJs{V) 

s =1 
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where j7fe(V)’s are determined by the BVP ([2.ID and (12.2j) for fixed Ck s and IZk s. 
Of course, the relations of individual fluxes Jk s to V contain more information ([21]) 
but it is much harder to experimentally measure the individual fluxes J7fc’s. 

2.1.1 Electroneutrality boundary conditions 

In relation to typical experimental designs, the positions X = ao and X = b o are 
located in the baths separated by the channel. They are the locations of two electrodes 
that are applied to control or drive the ionic flow through the ion channel. Ideally, 
the experimental designs should not affect the intrinsic ionic flow properties so one 
would like to design the boundary conditions to meet the so-called electroneutrality 

n n 

y z s c s = o = y z s n s 

S=1 S=1 

The reason for this is that, otherwise, there will be sharp boundary layers which cause 
significant changes (large gradients) of the electric potential and concentrations near 
the boundaries so that a measurement of these values has non-trivial uncertainties. 
One clever design to remedy this potential problem is the “four-electrode-design”: 
two “outer-electrodes” in the baths far away from the ends of the ion channel to 
provide the driving force and two “inner-electrodes” in the baths near the ends of the 
ion channel to measure the electrical potential and the concentrations as the “real” 
boundary conditions for the ionic flow. At the inner electrodes locations, the elec- 
troneutrality conditions are reasonably satisfied, and hence, the electric potential and 
concentrations vary slowly and a measurement of these values would be robust. We 
point out that the geometric singular perturbation framework for PNP type models 
developed in [T2[ [131 [22, 25) [26l [27] can treat the case without electroneutrality as¬ 
sumption; in fact, the boundary layers can be determined by the boundary conditions 
directly. 

2.1.2 Electrochemical potentials 

The electrochemical potentials Hk{X) consists of the ideal component and the 
excess component //“ where the ideal component 

/4 d (X) = z k e 0 <f>(X) + k B T In (2.4) 

is the point-charge contribution where Co is a characteristic concentration, and the 
excess component /x^( x ) accounts for ion size effects. As explained above, although 
not totally physical for ion channel problems in general, we will consider only the 
ideal component in this work that can act as guidance for further studies of more 
accurate models with excess component. 

2.1.3 Permanent charges and channel geometry 

The permanent charge Q{X) is a simplified mathematical model for ion channel 
(protein) structure. It is determined by the spatial distribution of amino acids in the 
channel wall, the acid (negative) and base (positive) side chains, more than anything 


(2.3) 
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else (IS!)- We will assume Q(X) is known and take an oversimplified description to 
capture some essence of its effects. For this paper, we take it to be as in (P2J), for 
some ao < A < B < bo, 


0i X £ [a 0 , A) L) (B,b 0 ] 
Q[X) 1 2Q 0 , Xe(A,B). 


(2.5) 


We will be interested in the case where |Qo| is large relative to /Vs and TVs. 

The cross-sectional area A(X) typically has the property that A(X) is much 
smaller for X £ ( A,B ) (the neck region) than that for X ^ ( A,B ). It is interesting 
to note that, in [23| . the authors showed that the neck of the channel should be 
“narrow” (small A(X) for X £ (A, B)) and “short” (small B — A) to optimize an 
effect of permanent charge. 


2.1.4 Dielectric coefficient and diffusion coefficients 

We assume that 


e(X) = e r is a constant, and T>k(X) = V(X)T>k (2-6) 

for some dimensionless function T>(X) (same for all k) and dimensional constant TV 
Note that the assumption ’D^(X) = T>(X)T>k is equivalent to the statement that 
T>k{X)/Vj{X) is a constant for k / j. Roughly speaking, the assumption says that, as 
the environment varies from location to location, its influences on the two diffusion 
coefficients T>k(X) and T>j(X) at the same location x are the same; that is, the 
two diffusion coefficients vary from one common environment to another common 
environment in a way so that their ratio is independent of locations. This is not a 
justification of this assumption but only an explanation of what it reflects. 

2.1.5 Main assumptions 

In the sequel, we assume the boundary electroneutrality condition in \2. A) , the form of 
the permanent charge in constant dielectric coefficient and diffusion coefficient 

property in 12. 61) . and the electrochemical potential is ideal pk = M/if in There 

are two key assumptions for this work to be discussed in terms of dimensionless 
variables below. 


2.2 Dimensionless form of the quasi-one-dimensional PNP model 

The following rescaling (see m) or its variations have been widely used for conve¬ 
nience of mathematical analysis. 

Let Co be a characteristic concentration of the ion solution. We now make a 
dimensionless re-scaling of the variables in system (12.111 as follows. 


2 _ ErEoksT 

eo(b 0 - a 0 ) 2 Co 
D(x) = V{X), 4>{x) 


, x = ——, h(x) 

bo — ao 

= <*(s) 


A{X) 

(bo - ao) 2 ’ 
Ck(X) T 


Q(x) 


Co 

Jk 


(bo — ao)CoT>k 


(2.7) 
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The dimensionless quantity Q(x) from the permanent charge Q in (12.51) becomes 


( 0, x € [0, a) U (6,1] 

\ 2Qp, x € (a, b), 


( 2 . 8 ) 


where 


0 < a = 


A- a 0 

bp — ap 


< b 


B -dp 
bp — dp 


< 1 , 


and, in terms of the dimensionless quantities, the subinterval (a, b) C (0,1) corre¬ 
sponds to the neck region [A, B\. 

Key assumptions. The parameter e is small and the parameter Qp > 0 is large. 

The case where Qp < 0 with |Qo| large can be treated in the same way. The 
assumption of smallness of the dimensionless parameter e is reasonable particularly 
because we are interested in the limit of large |<5o|- Thus, we can take large but 
fixed characteristic concentration Cp in the definition of e (12.71) . For example, if 
bp — dp = 25 nm and Cp = 10M, then the dimensionless parameter e ~ 10' - 3 m)- 
The mathematical consequence of smallness of e is that the BVP (12.91) and (12.101) 
can be treated as a singularly perturbed problem. A general geometric framework for 
analyzing the singularly perturbed BVP of PNP type systems has been developed in 
[EJ [Ml E3 So] for classical PNP systems and in [22] [25] [29] for PNP systems with 
finite ion sizes. 


In terms of the new variables in (12.71) . the BVP (12.11) and (12.21) becomes 


d 


h(x) dx 

dJk 
dx 


h(x 


dcj) 


dx 

= 0 , — Jk = 


= ~^2z s c s - Q(x), 


s=1 




with boundary conditions at x = 0 and x = 1 

0(0) =V, Cfc(0) = L k \ 0(1) = 0, Cfc(l) = R k , 


(2.9) 


( 2 . 10 ) 


where 


V := 


ep 

k B T 


V, 



R k ■■= 


Rk 

Cb‘ 


In this work, we will consider the BVP (12.91) and (12.101) for ionic mixtures with 
one cation of valence z\ = 1 and an anion of valence z% = —1. We will be interested 
in properties of ionic flow for large |Qo|- It turns out 


a = 


H{a) 


and f3 = 


m 


where H{x) = 


/' 

J o 


1 


-ds 


H( 1) ^ H( 1) v ' J o D(s)h(s) 

are the key parameters that characterize the effect of channel geometry and location 
of permanent charge. The physical meanings of these parameters could be seen from 
the special case when h{x) = hp and D{x) = Dp are constants. In this case, 


H(x) = 


Dphp 


which is proportional to the (scaled) length x of the region over [0, x] of the channel 
(conductive material), inversely proportional to the (scaled) cross-sectional area hp of 
the channel and to the (scaled) diffusion coefficient or electric conductivity Dp ; that 
is, H(x) is the resistance of the portion of the channel over [0, x\. 
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3 Approximations of fluxes [45J. 


We now recall the results on approximations of (</>, ci, C2, Ji, J 2 ) and nfl s from [45 J 
for the case where z\ = 1 and z 2 = — 1 with L\ = L 2 = L and R\ = R 2 = R. 

Based on the assumptions that e > 0 is SMALL and Q 0 > 0 is LARGE, one has 
expansions of fluxes in e and in v = l/Qo near e = 0 and v = 0. 

First, one expands the fluxes in e as 


Ji(e; v) = Ji{v) + O(s) and J 2 (s) = J 2 {iz) + 0(e), 

where Jfc(i'), depending also on the parameters (' V,L,R,H(l),a,(3 ), are the zeroth 
order in e terms of the fluxes. Then, one expands Jk{n) about v as 

J\{y) = J w + July + 0(v 2 ) and J 2 (v) = J20 + Jnv + 0{v 2 ). 


Thus, 

Ji(e; v) = J 10 + Jut' + 0(u 2 ,e) and J 2 (e; v) = J 20 + J 21 V + 0(u 2 , e). (3.1) 

Note that J 10 and J 2 o, depending on (V, L, R, H(l),a, ft), contain the leading order 
effect of the small v (or equivalently, large Qq). 

The following result is established in [ 15] . 

Proposition 3.1. One has 


Ji 0 — 0 , 

•h 1 = 


(1-/3 )L + aR 


{e V L - R); 


J20 — 


■hi = 


-a) - p)VeVL + aVR. 

2VLR 1 


H(l) (1 - P)y/L + aVIFVR 

(1 — j3 )L + aR 


(V e~ v L — VR), 


2H(l)(fl - a){Ve~ v L - y/R)(( 1 - /3 )VL + aVe~ v R) 3 
{—2(/3 - a) 2 (Ve~ v L - VH) 2 Ve~^LR 

+((1-/3 )L + aR){L-Re~ v ) (1 ~ ^ + aR VF^ln -±- 

-((1 - P)VL + aVe~ v R){Ve~ v L - VR) }. 


(3.2) 


A distinct implication of J 10 = 0 in Proposition 13.11 is that large (positive) per¬ 
manent charge Qo (or small v = 1/Qo) inhibits the cation flux. We will provide a 
detailed discussion in Section [4] on what happens to the internal dynamics that is 
consistent with this conclusion. 

To this end, we recall another immediate consequence of (13.2D (see [35] for more). 


Corollary 3.2. [Current Saturation] For large permanent charge Q 0 (small v = 
\/Qo) and to the leading order terms in v, each individual fluxes Jk’s, and hence, the 























total current I saturate as \V\ —>• oo; more precisely, one has 


lim J 20 

V —>+00 


lim J\i 

V —>+00 


lim Jn 

V — ¥—OO 


1 2 R 


lim J 2 i = 

V —>+00 

lim J 2 i = 
V— 7 —00 


I' 7 1 2L 

lim J 2 o = -777TT) 

1 ((1 - p)L + aR) 2 

(l-/3)2 2tf(l)(/?-a) 

1 ((1 - p)L + aR) 2 
~a 1 2H{l)(P-a) 


Proof. All the above (finite) limits can be derived from (|3.2I) easily. 
Note this is NOT the case if the permanent charge is small ([23]). 


(3.3) 


□ 


4 Internal dynamics for J 10 = 0. 

It follows from the Nernst-Planck equation in (12.91) that 

Jl r D l{X )mc, {X ) dx = '‘ l(0) - w(1) ' 

Thus, J\ has the same sign as that of /ii(0) — ^i(l); in particular, if /ii(0) — //i(l) / 0, 
then J\ / 0. 

For the setup of this paper, there are three regions of permanent charge Q(x ): 
Q(x) = 0 for x € [0, a) and x € ( b , 1] and Q(x) = 2Qq for x € [a, b] with large Q 0 or 
small v with v = 1 /Qq. A major consequence of large Q 0 is, from Proposition 13.11 
up to the leading order in v = 1/Qo an d in e, the flux of cation is J 10 = 0, even if the 
transmembrane potential 0. We will reveal the internal dynamics that 

lead to this conclusion. To do so, we will discuss what happens over each subinterval 
based on the approximated (of zeroth order in e) functions of profiles. Let 

(0(x; e, v),c k (x; e, v), J k (e, v)) = (0(x; v),c k (x; u), J k (u)) + O(e) (4.1) 

be the solution of the BVP (12.91) and (|2.10|) . For v > 0 small, one has the following 
expansions 


Ji(z/) =Jio + Juv + 0(v 2 ), 

where J 10 = 0 and J\\ are given in (13.21) . We will also provide figures for the profiles of 
rescaled electrical potential <f(x] is), cation concentration c\(x] is) and electrochemical 
potential n\(x\is) of cation, respectively. The parameter values used are 

e 0 =1.6022 x 10" 19 (C), k B = 1.3806 x 10' 23 (JA'" 1 ), T = 273.16 (K), 

V =0.01 (V), C = 0.5 (Af), 11 = 10" 5 (M), Q 0 = 10 8 (M), C 0 = 1 (M), (4.2) 

ao =0, bo = 10 (nm), A = 1/3 (nm), B = 1/2 (nm). 

In terms of the dimensionless parameters, 

V « 0.425, L = 0.5, R = 10~ 5 w 0, a = \,b = \, u = 10" 8 w 0, 

O Zj 
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The scaled cross-sectional area h(x) is taken to be 


{ 7 r(-x + r 0 + a) 2 , x€(0,a) 

7T7-q, x € (a, b ) 

^ t(x + r 0 - b ) 2 , a: € (6,1). 

The same setup will be used for figures in Section [5j 

Notice that the whole interval is (0,1), and it is divided into three subintervals 
(0, a), (a, b) and (6,1). Permanent charge is located over (a, b). The radius of the neck 
of ion channel is c. In the figures, we let ro = 0.5. We make plots of each function in 
interval (0,1), also in every subintervals. 

4.1 Internal dynamics over the interval (0,a) 

The leading order terms of (4>,c\) in (14. ip are derived in [35]. One has 
Proposition 4.1. For x € (0, a), 

4>{x] u) =<J)o(x) + 4>i{x)v + 0(v 2 ) with 

J20 TT ,^\ j. /_a C11O) , 2 J n ln ci 0 (x)_ 


Mx) =V - In (l (j>i(x) = - Cn * + 

V 2L / ciox) 


cio (a:) J20 L 

cio(x) =L - ^j-H(x), cn(x) = - Ju + J21 H(x)\ 

/iio(x) =/rio(0) = V + lnL, Hn{x) = 2^- In Cl0 j x ^ . 

020 L 

Figure 1 shows the profiles of ci(x; v), 4>(x; v) and ^i(x; u) over the interval (0, a). 
The following is then a direct consequence. 

Corollary 4.2. Over the interval (0,a), cio(x) = 0(1) and 3lo( x ) = 0> anc ^ hence, 
one has J\o = 0. 


4.2 Internal dynamics over the interval (a, b) 

It follows again from [45] that one has the following approximations. 

Proposition 4.3. For x € (a, 6), 

4>(x] v) = — In v + (j) o(x) + 4>i(x)v + 0(u 2 ) with 
4>o{x) = lnR — 21nl?o + In 2, 

4>i(x) =cf>l - A 0 + -^L(H(x) - H(a )); 
cio(x) =0, c n (x) = - J u (H(x) - H(a)), 

/Uio(x) = In R - 2 In B 0 + In 2 + In ^ 7^0 ~ -hi ( H{x ) — FI (a)) ), 
Mil (x) =(f>i ~A 0 + -^-(H(x) - H(a)), 
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Figure 1: Profiles of ci(x\u), 4>(x\i/), and n\(x\v) over interval [0, a]. Note that 
!j!\q(x) = 0 but fi\ (x; v) = ji\ \(x)i> + 0(v 2 ) / 0 as shown in the figure. 


where 


A 0 — 
B 0 = 
A\ = 


B i = 


ri = 


Ve^L 


(1 - + aVR 

Vr 

(1 - /?)v / e F Z + aVR 


((1 -P)L + aR), 
((1 -P)L + aR), 


b q l 


2 (P ~ a) 2 (L - A 0 ) 2 - a 2 (A 2 - B g) In ^ 

4(/3 - a)((l - /3)L + ai?)(L - A 0 ) 


A 0 i?o, 


(! ~ 0) f 2 (/3 - a) 2 (L - A 0 ) 2 - a 2 (A 2 - B 2 ) In 


4a(/3 - a)((l - /3)L + ai?)(L - A 0 ) 


ln^ 

In Sak 
111 A 0 R 


2^ — — ~T~ 

\B 0 A 0 


+ 


(3 — a 


a 


B, 


(L — Ao)) — 2—-(L — Ao) + L. 

Bq a 


AnBi 


O-DQ) 


Figure 2 includes the profiles of ci(x; v), 4>(x; v), and fii(x; v) over interval [a, b\. 
One has the following immediate consequence. 


Corollary 4.4. Over the interval (a,b), cio(x) = 0, and hence, J\q = 0. Further¬ 
more, /iio(a) = V + In L and fiio(b) = Ini?. As R -A 0, /^io(a) — /iio (b) -A oo. 

Proof. It follows from Proposition 14.31 directly. □ 
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Figure 2: Profiles of ci(x; v), <f>(x; u), /ii(x; v) over interval [a, 6], Note that cio(x) = 0 
but ci(x; v) = cio(a:) + cn (x)v + 0(v 2 ) = c\\{x)u + 0(y 2 ) / 0 as shown in the figure. 
Also note the large drop of ^i(x) over this interval due to that R is taken to be small. 
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4.3 Internal dynamics over the interval (6,1) 
Proposition 4.5. For x € (6,1), 


4>(x] v) =cf>o(x) + 4>i(x)is + 0{v 2 ) with 
4>o(x) = In R — In cio (x), 


fr(x) =4 - I ( B $- - gocn(^) - fiCioM + 2Jn ln£ io(£) ; 

2 \ Bo J Bqcio(x) J20 Bo 

c w (x) =Bo - ^(H(x) - H(b)), cn(x) = B 1 - — + hl (H(x) - H(b )); 


Vio(x) =Mio(l) = In R, 


/in(x) = (x) + 


cn(x) 

cio(x)’ 


where 


b ln ^ 

4)1 

lU A^R 


- — 

\Bo A 0 


+ L^ {L _ Ag) \ 2 I ?A + Bo. 

a J B 0 


Figure 3 shows the profiles of ci(x; v), <p(x] v), and ^i(x; v) over interval [b, 1]. 





Figure 3: Profiles of ci(x;i'), and gb\{x\v) on interval [6,1]. Note that 

// 10 (x) = 0 but ^(x-u) = /iio(x) + li' u (x)v + 0(v 2 ) = /i' n (x)i/ + 0(is 2 ) / 0. 

Corollary 4.6. Over the interval (6,1), noting Bq = 0{VR), cio(x) = 0(1) and 
Li[(x) = 0, and hence, J 10 = 0. 
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4.4 Summary of mechanism for ,J\q = 0. 

From the above discussion, we conclude that the mechanisms for J\q = 0 are different 
over each subintervals. More precisely, one has that, 

• over the first interval (0,a), Jio = 0 is the result of constant electrochemical 
potential fi\o(x) = /iio(0) = ^io(a) (so that // 10 (x) = 0) while cio(x) ^ 0; 

• over the last interval (6,1), similar to that over the first interval (0, a), Hio(x) = 

/iio(l) so that = 0 while cio(a:) / 0, and, for R small, c\o(x) = 0(\TR) 

over this interval; 

• over the middle interval (a, b ) where permanent charge is not zero and large, the 
electrochemical potential ^io(x) is not a constant so that fi' 10 (x) / 0, however, 
cio(x) = 0, in particular, the drop of ^io(x) over this interval equals the trans¬ 
membrane electrochemical potentials, that is, /Jio(o) — /iio (b) = /ii(0) — ^ti(l). 

Here we provide the profiles of concentration (Fig. 4), electrical potential (Fig. 
5) and electrochemical potential (Fig. 6) over the whole interval [0,1]. 

The figures of ci(x) and cj)(x ) over interval (0, xo) are not continue, because we 
make the plots of system (12.11) with e = 0. For the limiting system at x = a and 
x = b, there are two fast layers, ci(x) and 4>{x) changes very fast, but ^\{x) keeps 
the same value in fast layers. Recall that = /ii(0) — ^i(l) = ksTiV + lnL — lnii), 
if R = 10" 5 , /4 ~ 2.5902 x 10" 23 . 


Concentration of Cation in (0,1e-8) 


0.8 


0.6 


(r" 0.4 


0.2 


0 


-0.2 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


X 


xIO " 8 


Figure 4: Profile of ci(x‘,u) over [0,1]. Note that cio(x) = 0 for x G (a, 6) and, 
for x € (5,1), cio(x) = 0 (VR) is not zero but small for R = 10 5 chosen for the 
numerics. 
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Figure 5: Profile of <f>(x\ v) over [0,1] 



x10" B 


Figure 6: Profile of ni{x\v) over [0,1]. Note the large drop of /ii over the interval 
(a, b) = (1/3,1/2) and the nearly constant values over the other two subintervals. 
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5 Declining phenomenon and internal dynamics 


In this section, we will show that large permanent charge is a mechanism for the 
declining phenomenon described in the introduction. Recall that, by the declining 
phenomenon , we mean the following observed experimentally. 

For fixed V and L\ = L 2 = L, as R\ = R 2 = R decreases to zero, the flux of 
counterion (J 2 in the setting since Qo > Oj decreases monotonically to zero. 

This phenomenon is rather counterintuitive since the transmembrane electrochem¬ 
ical potential for the counter-ion 

/12(C)) — ^2(1) = zfiV + RL2 — lni?2 —> +00 as R2 = R —>• 0 . 

5.1 Experimental phenomena are consistent with our analysis 

The result in (13.21) actually justifies the observation, up to the leading order J 20 hi v 
for v near 0 (or for Q 0 —>• + 00 ); that is, 

Proposition 5.1. The leading order term J 20 of the flux is monotone and concave 
downward as a function of R and, as R —>• 0, J 20 —> 0. 

Proof. Indeed, from the expression of J 20 in (13.21) and treating J 20 as a function of w 
where R = w 2 for convenience, one has 

VL Ve~ v Lw — w 2 
277(1) (1 — fl)y/L + aVe~ v w 

It is clear that J 2 o(uj) —> 0 + as w —>• 0 + . Note that the derivative of J 20 in w is 

, 1 (1 — fl)LiV e~ v — 2(1 — j3)>/Lw — a\J e~ v w 2 

2H(1) [(1 — fl)VL + aV e~ v w ] 2 

Thus, from the expression of the numerator, 

if w is smaller than some wq, then J' 20 (w) > 0, 

and hence, as w —>• 0 + (or equivalently, 1?2 = R —>• 0) monotonically over the interval 
[0, rco], J 2 o(w) —>• 0 monotonically. 

It is not hard to show that for R 2 = R > 0 but smaller than some positive value, 
the graph of J 20 as a function of R is concave downward. □ 

In Figure 7, the horizontal axis is for R and the vertical for J 20 . We fix V = 0.425 
and L = 0.5 as in ()4.2I) . and vary R £ (0,10~ 3 ]. The monotonicity and concave 
downward features of the graph are apparent. 

Remark 5.1. We comment that, for large permanent charge, the declining curve phe¬ 
nomenon occurs when the transmembrane electrochemical potential 2 ( 0 ) — / 12 (f) is 
increasing to infinity in a particular way; that is, as R —» 0. If one increases the trans¬ 
membrane electrochemical potential /^ 2 (0) — / 12 (f) in a different manner, for example, 
as \V\ —>■ 00 or as L —>• 00 , Corollary 13.21 shows that the declining curve phenomenon 
does not happen. 


f6 












R xio' 2 


Figure 7: Declining curve: J 20 vs R for R € (0,10 3 ] with L = 0.5 and V = 0.425. 

It is also important to note that, when the next order term J 2 \v is considered, 
then, as R 2 = R — >• 0, J 21 — >• 00 . But, if R — >• 0 and u In R, — >• 0, then J 2 \v — >• 0. 

Thus, only when v is very small (Q 0 is very large), is the term J 2 \v not significant, 
and hence, the term J 20 dominates the described behavior. □ 


5.2 Mechanism of declining phenomena from the profiles 

Recall, from the Nernst-Planck equation in (|2.9I) that 

~’h = -j^D 2 (x)h(x)c 2 (x)-^ia 2 (x). 

Since D 2 (x) and h(x) are fixed, we will treat them as of order 0(1) quantities so that 
they do not contribute much to the near zero flux scenario as R —>• 0. Thus, as far as 
the near zero flux mechanism is concerned, one has 


-J 2 


C 2 (x)—/x 2 (x). 


(5.1) 


One sees that the gradient of the electrochemical potential i 2 (x ) is the main 
driving force for the flux J 2 . Intuitively, large drop of (or transmembrane) electro¬ 
chemical potential /a 2 (0) — fa 2 (l) of fi 2 produces large flux J 2 . In this sense, the 
declining curve phenomenon is rather counterintuitive. A careful look at (15.11) reveals 
that there is only one possibility for the declining curve phenomenon; that is, when¬ 
ever -^/i 2 (x) is large, c 2 (x) has to be much smaller in order to produce a small flux 
| J 2 1. We will apply the analytical results of the internal dynamics from ([55]) to show 
that this is indeed the case. 

For our setup, there are three regions of permanent charge Q(x): Q(x) = 0 for 
x € [0, a) and x G ( b , 1] and Q(x) = 2Qq for x €E [a, b] with large Qo or small v with 
v = 1/Qo- 

For fixed V and L, fj, 2 (0) — /i 2 (l) = z 2 V + In L — In R « - Ini? > 1 for small R. 
We need to understand 


17 





(i) HOW the electrochemical potential /j , 2 drops an order 0 (— Ini?) 1 over the 
interval x € [ 0 , 1 ]; 

(ii) HOW J2 can be small, for small 12 and R, at every x € [ 0 . 1 ]; 

(iii) Most importantly, HOW the above two things, with the constraint ( 15 . 11 ) . can 
happen simultaneously. 

There are two small parameters 12 and R in our considerations. The relative sizes 
of these two parameters is relevant for the result. We will assume 12 In R -C 1 . 

We now discuss what happens over each subinterval based on the approximated 
(of zeroth order in e) functions of profiles. To do so, let 

(cj)(x- e, u),c k {x\ e, zz), Jk(e, 12)) = (0(2; u),Ck(x-, 12 ), + O(e) 

be the solution of the boundary value problem. For v > 0 small, one has the following 
expansions 

J2 (y) = J20 + J 2 \v + 0 (is 2 ), 

where J20 and J21 are given in ( 13 . 2 j) . The expansions in 12 for <f(x-, 12), c\ (2;^), and 
02(2; 12) are not regular and are qualitatively different over the subintervals ( 0 , a), 
(a, b) and ( 6 , 1 ). They will be given explicitly in each subsection below for us to 
understand what happens over each subinterval. 


5.2.1 Internal dynamics over the interval (0, a) 

The leading order terms of (4>, C 2 ) are derived in [25]. One has 
Proposition 5.2. For x € (0, a), 

4>(x; 12 ) =(^ 0 ( 2 ) + cfi(x)i 2 + 0(v 2 ) with 

/ / \ -rr i ( -1 J 20 rj/ \ \ , ( V C 2 l( 2 ) 2 J U C 20 (x) 

Mx) =V - In (1 - —H{ 2 )j , Mx) = + -j- In — ; 

C 20 {x) =L - -^-H(x), c 2 i (2) = - Jn + 2 J21 H{ 2); 

V2o(x) = -V + 2 In 020(2) - InL, fx 2 i(x) = 2 (-ihl _ 'Ill i n Q(> ^ l \ . 

\ C20 (x) J20 L ) 

Figure 8 shows profiles of c 2 (x\v), (f{x',v) and ^ 2 ( 2 ;^) over the interval (0,a). 
Corollary 5.3. Over the interval (0,a), 02 ( 2 ;^) = 0(1) BUT 


fj,o(x; v) = — . 20 . + 0 (n In R) = 0 (Vr , 12 In R). 

h\X)c 2 o (2) 


Therefore, from \ 5 . 1 \) . 


d 


J 2 « 02(2)—^2(2; 12) = 0(1)0(Vr, ulnR) = 0 (y/~R, 12 In R) 
dx 


and ^2(2) drops an order of 0 (y/~R, 12 In I?) over the interval ( 0 , 0). 

In particular, for zeroth order in 12, J 2 ~ J20 = 0 (y/~R) and p. 2 (x) = o{Vr ) over 
the interval ( 0 ,a). 
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Figure 8: Profiles of Ck(x',v), (j)(x',u), and Hk(x\ v) over interval [0, a]. 

5.2.2 Internal dynamics over the interval (a, 6) 

It follows from [45] that 
Proposition 5.4. For x € ( a,b), 

4>(x] u) = — In v + 4>o(x) + cf>i{x)v + 0(v 2 ) with 
(j)Q (x) = In R — 2 In Bq + In 2, 

M x ) =0? - A 0 + ^j-(H(x) - H(a )); 

c 2 (x;i/) - J n (H(x) - (a))^ v + 0(^ 2 ); 

h 2 o(x) = — In R + 2 In So — In 2, 

F2i{x) =-01 +Ao- ^-(H{x) - H(a)), 
h22(x) =^Al - Jn(H(x) - H(a )); 
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where 


Ao 

A\ 

Bi 




V / e r L((l - p)L + aR) R _ y/R((l - P)L + aR) 
(1 — /3)V e v L + ay/H (1 — (3)V e v L + ay/R 

2 (P - a ) 2 ( L ~ A 0 ) 2 ~ a 2 (Al - Bl) In 

4(/3-a)((l - P)L + aR)(L - A 0 ) 0 °’ 


(1 - P) (2(/3 - af(L - A o y - a 2 (A 2 - 5 2 ) In g$) 


lnffi 


4a(/3 - a)((l - 0)L + aR)(L - A 0 ) 

■(Si _ A 


A)-Bo, 


(Ki 1 - t) + — ( L - A))) - a|i - —(L - Ao) + L. 

\ \ B 0 A 0 / a J B 0 a 


The profiles of C 2 (x; v), <f>(x\ v), and ^(x] v) are shown in Figure 9. 




Figure 9: Profiles of C 2 (x;i'), (ffais), and H 2 (x-,v) over interval [a, b] 


Note that, over this interval, C 2 (x; v) is singular in v. To understand what happens 
to the internal dynamics for anions over the interval (a, b), we need to resolve this 
singularity by considering at least 0(i/)-terms. 

Corollary 5.5. Over the interval (a, 6), C 2 (x;u) = 0(l/v) 3> 1 and fi 2 (x',v) = 
0{lnR) > 1 BUT 


/4( x '; v) ~ /4i ( x ) u = = o(vVr). 

zh{x) 

Therefore, from 1,5. 1\) , 

J 2 k, 0(l/v)0(v'/R) = 0(y/R) and /i 2 (x;i/) only drops 0(v'fR) over ( a,b ). 
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Remark 5.2. Note that p 2 (x) drops much less over the interval (a,b) than its drop 
over the interval (0, a). But both drops are small and contribute nearly nothing to 
the total drop p 2 (0) — /^(l) = 0(— In R ) 1. The only way to realize the large total 

drop p 2 (0) — h 2 (l) is that /i 2 (x) drops A LOT over the subinterval (6,1). Indeed, this 
is the case as claimed below. □ 


5.2.3 Internal dynamics over the interval (6,1) 
Proposition 5.6. For x € (6, 1), 


4>(x\ v) =<fio(x) + (j>i(x)u + 0(y 2 ) with 
4>o(x) = In R - In c 2 o (x ), 


J. j,b 1 ( 0 2 2B \ - B 'o\ B 0 C 2 l{x) - B lC20 (x) 2Jn,_ C 2 0 (s). 

Mx) ^ “ 2 1 0 ■ —^) - FFF) - + IF 1,1 

C 20 (x) —Bq - F(H(x, - H(b)), C 2 i(x) - B J " ~ J '“ (fflx) - H(b)); 


[J ,20 (x) =2 In c 20 (x) - In R, fi 2 i{x) 


C 21 (x) 
C 2 o(x) 


M x ); 


where 


ln 7T 

In TBL. 
111 A 0 R 


\B 0 Aq 


+ ~——(L ~ A))'] _ + Ho¬ 
ck J Bq 


The profiles of C 2 (x; u), cj)(x] v), and [i 2 (x; v) over (6,1) are shown in Figure 10. 

Corollary 5.7. Over the interval (6,1), noting Bq = 0{y/R), c 2 o(a;) changes from 
C 20 (6) = Bq = 0(y/R) to C 20 (1) = R monotonically, and fj, 2 (x) changes from [ 12 (b) = 
0(1) to ^ 2 ( 1 ) = Ini?. 

Therefore, for x € (6, 1), from 15.il) . 


J 2 « 0{VR) and /r 2 (x) drops 0(— Ini?) over (6,1). 


Remark 5.3. Note that, over this interval, the order of [i' 2 {x ) varies in x from /z 2 (&) = 
0(1) to /u 2 (l) = 0(l/y/R) but overall drops is 0(— Ini?). This is different from what 
happened over the intervals (0, a) and (a, 6). □ 


5.3 Summary of mechanism for declining phenomenon. 

In summary, with the technical assumption that zzlnf? < y/R, we have, as i? -> 
0, J 2 = 0(y/R) over the whole interval (0,1) but with completely DIFFERENT 
scenarios over different subintervals (0, a), (a, 6) and (6,1). More precisely, 

(i) over the subinterval (0,a), one has c 2 (x) = 0(1) but /i 2 (x) = 0(y/~R) <C 1 
so that J 2 = 0(y/R)\ (Note that the drop of /r 2 over the interval (0,a) is of 
order 0 (V~R), which has nearly no contribution to the drop of p 2 over the whole 
interval (0,1).) 
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Figure 10: Profiles of c 2 (x;v), (j)(x',u), and p 2 (x;is) over interval [6,1] 


(ii) over the subinterval (a, 6), C 2 (x) = 0(l/v) but p' 2 {x) = 0(vyfR) < 1 so that 
J 2 = 0(VR)-, (Note that the drop of p 2 over the interval (a,b) is of order 
0(uy/R), which is even smaller than that over the subinterval (0, a) and, of 
course, has nearly no contribution to the drop of p 2 over the whole interval 

(o,i)0 

(iii) over the subinterval (6,1), different from what happened over each of the previ¬ 
ous two subintervals, the orders of C 2 (x) and p 2 (a :) are NOT uniform in x £ (6,1) 
but the drop of p 2 (x) of 0(ln R) is fully realized over this subinterval (see Re¬ 
mark 15.31) . 

Here we provide the profiles of concentration (Fig. 11), electric field (Fig. 12) and 
electrochemical potential (Fig. 13) of the anion over the whole interval [0,1]. The 
figures of Ck(x) and <f{x) over interval (0,xo) are not continue, because we make the 
plots of system {2TD with e = 0. For the limiting system at x = a and x = b, there 
are two fast layers, q,(x) and 4>(x) changes very fast, but Pk( x ) keeps the same value 
in fast layers. Recall that p 2 = ^ 2 ( 0 ) — ^ 2 ( 1 ) = ksT^—V + InL — InR), if R = 10 -5 , 
then P 2 ~ 2.4780 x 10~ 31 , that’s why in the last figure, P 2 ( x ) can n °t reach the 
infinity. 
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Figure 11: Profile C 2 (x;u) over [0,1] 



x ■; 10" B 


Figure 12: Profile of (j)(x;u) over [0,1] 



xIO ' 8 


Figure 13: Profile of H 2 {x\u) over [0,1] 
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6 Concluding remarks 


In this work, we examine effects of large permanent charges on ionic flow through ion 
channels based on a quasi-one dimensional Poisson-Nernst-Planck model. We show 
that one of the defining properties of transporters, obligatory exchange, can arise in 
an open channel with just one structure. When the permanent charge is large, the 
current carried by counter ions, majority charge carriers with the opposite sign as 
the permanent charge, can decline, even to zero, as the driving force (the gradient 
of electrochemical potential) increases. We also show that large permanent charges 
essentially inhibit the flux of co-ions, regardless of the magnitude of transmembrane 
electrochemical potential. 
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